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Data from the CCFR E770 Neutrino Deep Inelastic Scattering (DIS) experiment at Fermilab contain events 
with large Bjorken x (x > 0.7) and high momentum transfer (Q 2 > 50 (GeV/c) 2 ). A comparison of the 
data with a model based on no nuclear effects at large x, shows a significant excess of events in the data. 
Addition of Fermi gas motion of the nucleons in the nucleus to the model does not explain the excess. 
Adding a higher momentum tail due to the formation of "quasi-deuterons" makes some improvement. 
An exponentially falling F2 oc e~ sfx ~ x °^ at large x, predicted by "multi-quark clusters" and "few-nucleon 
correlations", can describe the data. A value of s = 8.3 ± 0.7(stat.) ± 0.7(sys.) yields the best agreement 
with the data. 
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Deep inelastic scattering (DIS) neutrino interactions 
are ideally suited for measuring the nuclear structure 
functions over a wide range of kinematic variables. The 
CCFR experiment has collected over 10 6 charged current 
events where neutrinos and antineutrinos scattered off of 
an iron target. We have recently published yj] structure 
functions for the Bjorken x region of x < 0.75 based on 
these data. In this paper, we extend the measurement 
into the x > 0.75 region We have approximately 
2000 events with x > 0.75. 

In the infinite momentum frame, x is the fraction of 
the nucleon's momentum carried by the struck quark. As 
such, x is kinematically constrained to be less than 1 . In 
the same model, quark counting rules || would imply 
that the quark structure functions should decrease like 
(1 — x) 3 near x = 1. 

When x is measured in DIS experiments, it is assumed 
that the incident lepton scatters off a nucleon at rest in 
the laboratory. Under this assumption and neglecting 
terms proportional to the nucleon's mass, x takes on the 
form 



MfjEhad 




where Ei is the energy of the incoming lepton, Ey is the 
energy of the outgoing lepton, is the laboratory angle 
between the two leptons, E^ad is the energy of the outgo- 
ing hadron shower, and Mm is the mass of the struck nu- 
cleon. However, for nuclear targets (like the iron detector 
of CCFR) the nucleons need not be at rest in the labora- 
tory. The nucleons themselves have Fermi motion. Fur- 
thermore, the quarks inside a given nucleon may be af- 
fected by the surrounding nucleons either through the ex- 
change of mesons inside the nucleus, through the forma- 
tion of few-nucleon states, or through the condensation 
of the nucleons into multi-quark clusters ( "bags" ) . All of 
these phenomena have been used to explain the EMC ef- 
fect Q. Nuclear motion and quark interactions "smear" 
out the nucleon's structure functions, taking events from 
one region in x and Q 2 and moving them into adjacent 
regions. In the regions below x = 0.75, the differences 
between the nuclear and nucleon structure functions are 
small mainly because the nucleon structure functions are 
slowly varying. At x = 0.7, the nuclear correction is ex- 
pected to be about 10% |T|. However, in the higher x 
region, nucleon structure functions must go to zero while 
the nuclear structure functions need not. In fact, most of 
the cross section near x = 1 and all of the cross section 
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for x > 1 comes from the fact that the quark exists in the 
nuclear environment. Thus, the high x region is ideally 
suited for investigating nuclear effects in DIS scattering. 

Event selection in this analysis is very similar to previ- 
ous CCFR analyses Jl]||. However, because there are 
so few true high x events, we are very susceptible to 
background events from the lower x region being mis- 
measured and appearing as high x events. To minimize 
this effect, we further restrict the kinematic region for 
this analysis over previous CCFR analyses to the region 
where x is well measured. The resolution of x is poor 
at large muon energies, at small hadron energies, and at 
small muon angles. We also require that the hadronic en- 
ergy is well contained in the calorimeter and has a small 
measurement error and that the outgoing muon traverses 
the whole length of the toroid spectrometer. For these 
reasons, we require events to be in the kinematic region 
15 GeV < Ef, < 360 GeV, 20 < E had < 360 GeV, 
Q 2 < 400 (GeV/c) 2 , and M > 17 mrad. To mini- 
mize the uncertainties from higher twist contributions 
and from resonant scattering, we accept events where 
Q 2 > 50 (GeV/c) 2 . The average Q 2 of the accepted 
events is 120 (GeV/c) 2 . 

In addition to these quality cuts, we made various 
checks on the muon momentum and hadronic energy 
measurements. We compared different measurement al- 
gorithms and hand-scanned all of the high x events. 
There were about 30 events with an apparent discrep- 
ancy between two different hadron energy or muon en- 
ergy measurement algorithms, indicating possible mis- 
measurements. These events were removed from the sam- 
ple. 



x r 

c 

go.25 - 
W 



0.2 - 



0.15 - 



0.1 - 
0.05 - 

r " 

0.6 0.7 0.8 0.9 1 1.1 1.2 

Measured x 



FIG. 1. Error in a; as a function of measured x in the data. 



As with any steeply falling function, the measured 
distribution is distorted by measurement resolution. In 
CCFR, the hadron and muon energy resolution functions 
are measured from test beam data and from internal cal- 
ibrations They are well known over four orders of 
magnitude. An extensive Monte Carlo program has been 
written to simulate the effects of the detector on neutrino 
interactions. With this program, we predict the observed 
x distribution of a given input model. For models with 
free parameters, we vary the free parameter until we get 
the best agreement between the predicted x distribution 
and that observed. The x resolution as a function of x 
is shown in Fig. |l|. The number plotted is the resolution 
averaged over the accepted Monte Carlo events. 

Neutrino DIS can be completely parameterized (up to 
effects that are proportional to the muon mass over the 
neutrino energy) by three structure functions: 2xF\, F%, 
and xF%. In the quark parton model, 2xF\ and F2 dif- 
fer by the distribution of effective scalar partons in the 
nucleon which arise from the apparent transverse mo- 
mentum of the quarks. At high Q 2 , the scalar effects 
are expected to be insignificant. 2xFi and xF^ differ 
by a contribution from sea quarks. At high x 7 the sea 
quark contribution to the structure functions is also neg- 
ligible. Therefore, for this analysis, we assume that [||, 
for x > 0.6, 

2xF 1 (x, Q 2 ) = F 2 (x, Q 2 ) = xF 3 (x, Q 2 ). (2) 
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FIG. 2. (a) Measured x distribution. Comparison of 
measured x distribution and (b) distribution predicted by 
Buras-Gaemers structure functions and by CTEQ4M struc- 
ture functions, (c) distributions predicted by a flat Fermi gas 
model and by the Bodek-Ritchie model, (d) distributions with 
an exponentially falling F2 with s=8.3. All error bars repre- 
sent only the statistical errors. 

As stated above, simple nucleon models of neutrino- 
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iron scattering do not include the nuclear environment 
and therefore are not expected to describe our x distri- 
butions. However, it is instructive to compare a model 
including only nucleon structure functions with our data 
to see the difference. Fig. ^(a) shows the measured x 
distribution from the data. In Fig. ||(b) , we show our ex- 
pected x distribution generated from the Buras-Gaemers 
(BG) and from the CTEQ4M parameterizations of the 
structure functions The former have been fit to our 
x < 0.7 data and, as such, include the nuclear effects that 
are seen at lower x. Likewise, the latter uses CCFR data 
in their fits and therefore, to some extent, must include 
nuclear effects. However, both parameterizations are ap- 
proximately proportional to (1 — x) 3 as x goes to 1 and 
therefore do not include high x nuclear effects. As can 
be seen in Fig. |(b), the BG and CTEQ4M parameteri- 
zations seriously underestimate the high x cross section. 

When a nuclear Fermi gas model with a flat momentum 
distribution up to a Fermi surface of 257 MeV/c is added 
to the BG parameterization, the result again underesti- 
mates the high x measurement as shown in Fig. ||(c). It 
has been found in low energy electron-nucleus scattering 
measurements and in analyzing the EMC effect in high 
energy experiments that a high momentum tail has to be 
added to the pure Fermi gas distribution to reproduce 
the data. The model of Bodek and Richie || introduces 
the high momentum tail through the inclusion of quasi- 
deuteron scattering. Fig. ^(c) compares this model with 
our data. While the comparison with our data is still not 
perfect, it is slightly better than that for the models with 
no nucleon motion or with only Fermi motion. 

After the discovery of the EMC effect, various more ex- 
otic nuclear models were proposed. One such model is the 
few-nucleon correlation model of Frankfurt and Strikman 
[ 10 1 . Other authors such as Kondratyuk and Shmatikov 
[11 1 suggest models where the neutrinos scattered off 
higher quark-count states (6-quark, 9-quark, etc.). In 
such models the individual quarks can have a higher ef- 
fective momentum than in the isolated nucleon. Both 
the few-nucleon correlation models and the multi-quark 
cluster models can be parameterized as an exponentially 
falling i<2 structure function (i^Oz) oc e~ sx ) in the large 
x region. We have used the parameterized structure func- 
tion from the BCDMS collaboration |l^] and varied the 
exponential slope parameter s to minimize the x 2 differ- 
ence between the Monte Carlo prediction and the mea- 
sured distributions in the region 0.6 < x < 1.2. We find 
that s — 8.3 minimizes the \ 2 at 8.2 for 12 degrees of 
freedom. The region that changes the total x 2 by 1 unit 
leads to an error estimate of 8s = ±0.7. Fig. ^(d) com- 
pares our best exponential fit to the data. Fig. |^ shows 
i*2 with its error band from this exponential model nor- 
malized to the x — 0.65 point from the CCFR structure 
functions [Rj. 
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FIG. 3. The exponentially falling F2 structure func- 
tion for Q 2 = 125(GeV/c) 2 along with the last two points 
from the previous CCFR analysis. The solid line corre- 
sponds to our best fit to an exponetial F2 structure func- 
tion. The area included by the dashed lines corresponds 
to the area ±lcr statistical error around the best fit. The 
dotted-dashed lines indicate the region allowed by adding 
both statistical and systematic errors in quadrature. The 
normalization of this graph is taken from all the data with 
0.6 < x < 0.7. The CCFR points plotted are only for the 
data 100 (GeV/c) 2 < Q 2 < 150 (GeV/c) 2 . 

We have investigated various systematic effects that 
might affect our determination of the exponential slope 
parameter. These are given in Table 1. The biggest ef- 
fect comes from our assumed muon energy resolution. If 
we broaden or narrow the resolution function from the 
measured one by one standard deviation, we can change 
the value of the slope parameter by ±0.6. Other system- 
atic effects, given in Table 1, lead to a total systematic 
error on the slope parameter of ±0.7. Thus our measured 
exponential slope parameter is 



s = 8.3 ± 0.7(stat.) ± 0.7(sys.). 



(3) 



Systematic Factor 


Error in s 


Energy scale 


0.2 


Relative calibration 


0.1 


Incoming neutrino angle 


0.1 


Outgoing muon angle 


0.1 


Hadronic energy mismeasurement 


0.1 


Radiative corrections 


0.1 


Muon resolution function 


0.6 



TOTAL SYSTEMATIC ERROR 



0.7 



TABLE I. Systematic errors. 
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Other experiments have measured the expontential 
slope parameter s for nuclear targets in the high x re- 
gion. Experiment E133 (u| at SLAC measured the low 
energy electron-aluminum DIS cross section. Their value 
of s = 7 — 8 cannot be directly compared with our result 
since their high x events come from the resonance region 
at low Q 2 and not from the DIS region. The muon-carbon 
data from BCDMS comes from the DIS region and should 
be comparable to ours. They find a value of s = 16.5±0.5 
(much steeper slope). This may indicate that the inter- 
nucleon effects are much greater for our iron target than 
for their carbon target, although this is not expected the- 
oretically H] . Recent studies at lower energies |u| from 
SLAC-NE3 and CEBAF E89-008 have suggested a con- 
sistent nuclear picture if £ = 2x/(l + + scal- 
ing is used. However, their exponent of approximately 
s = 17 at x = 1 is inconsistent with our data (in our Q 2 
regime, ^ « a;). Our result agrees with the theoretical 
prediction of s = 8 — 9 by Strikman and Frankfurt and 
is larger than the Baldin's prediction Jl(| of s = 6. 

The model that we used in this analysis assumed a 
Q 2 dependence of the structure function of ^2(2;, Q 2 ) oc 
(Q 2 ) -018 . This dependence is completely compatible 
with our data. 
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